Moving Observers, Non-orthogonal Boundaries, and 

Quasilocal Energy 



IS. Booth*and R.B. Mann* 
Department of Physics 
University of Waterloo 
Waterloo, Ontario 
N2L 3G1 

October 1, 1998 



Abstract 

The popular Hamilton-Jacobi method first proposed by Brown and York for defin- 
ing quasilocal quantities such as energy for spatially bound regions assumes that the 
timelike boundary is orthogonal to the foliation of the spacetime. Such a restriction 
is undesirable for both theoretical and computational reasons. We remove the orthog- 
onality assumption and show that it is more natural to focus on the foliation of the 
timelike boundary rather than the foliation of the entire four dimensional bound re- 
gion. Reference spacetimes which define additional terms in the action are discussed 
in detail. To demonstrate this new formulation, we calculate the quasilocal energies 
seen by observers who are moving with respect to a Schwarzschild black hole. 



1 Introduction 



Gravitational thermodynamics and its relationship to the Euclidean-action formulation of 
quantum gravity have been of increasing interest in recent years. This relationship was first 
explored by Gibbons and Hawking who argued that the Euclidean gravitational action is 
equal to the grand canonical free energy times the reciprocal of the temperature associated 
with a black hole (or cosmological) event horizon [0. A more recent extension of this work 
by Brown and York involved consideration of the formulation of the partition function for 
gravitating systems of finite spatial extent |J. Starting from a spacelike foliation of a 
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finite region of spacetime, and a timelike vector field defining a flow of time, they studied 
the Einstein-Hilbert action using a Hamilton- Jacobi type analysis. Decomposing the action 
according to the foliation and flow of time, they showed that natural candidates arose for 
quantities such as energy and momentum. These quantities were defined quasilocally, i.e. 
for a region of finite spatial extent containing a gravitating system. 

For a number of reasons this analysis and its associated quasilocal quantities have gen- 
erated much interest and found a multitude of uses. First, all physical systems with which 
we have any experience have a finite spatial boundary. Indeed one of the central concepts 
in thermodynamics is that of a system and a reservoir that are separated by a partition. 
Quasilocal quantities admit a physical realization of these concepts so that thermodynamics 
may be applied in a sensible way. As such, in the literature this analysis has been used ex- 
tensively in the study of black hole thermodynamics (for example in [|], |5|, |6|). Among other 
places, this work has found application in studies of the distribution of gravitational energy 
in a variety of spacetimes (for example 0) and also in examining the quantum mechanical 
creation of pairs of black holes (for example ||). A very similar Hamiltonian decomposition 
of the action has also been executed by Hawking and Horowitz || . 

However an acknowledged incompleteness exists in the quasilocal formulation in that 
(apart from two exceptions mentioned below) the spacetime foliation is always assumed to 
be orthogonal to the timelike boundary. While this is the case for many standard exam- 
ples (such as black holes surrounded by a set of stationary observers) it is a fairly strong 
restriction. For example, within the confines of this orthogonality assumption it is extremely 
difficult to calculate the quasilocal quantities seen by observers who are falling into a black 
hole. Furthermore when one considers variations of the metric (as one actually does dur- 
ing the quasilocal Hamilton- Jacobi analysis) the orthogonality assumption implies that the 
variations are not general, but instead restricted to those that preserve the orthogonality. 

The requirement that the timelike boundary of the finite region be orthogonal to its 
spacelike boundary was dropped by Hayward fLOf , who considered how the basic Hilbert 
action / should be modified so that solving 51 = for general variations of the metric (subject 
to the boundary condition that boundary metrics should be held constant) will produce the 
Einstein equations in the usual way. However this was from a purely Lagrangian viewpoint 
- no consideration was given as to how these variations would decompose in accordance with 



the spacetime foliation. That approach was recently considered by Hawking and Hunter [] 1 



and Lau [fL^| , who addressed the non-orthogonal situation from a Hamiltonian perspective. 

In ref. JTl| Hayward's action was broken down according to the foliation of the spacetime, 
a Hamiltonian proposed, and two sample calculations performed where the boundaries were 
non-orthogonal. However there was no attempt made in this treatment to consider the 
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variation of the action and show an agreement between the quasilocal quantities suggested 
by that approach and the direct Hamiltonian deconstruction of the action. Furthermore, in 
order to deal with the non-orthogonal intersections, the authors found it necessary to impose 
somewhat complicated restrictions on background comparison spacetimes. 



Lau's main interest in |L2| was to reformulate the quasilocal quantities of |J in terms of 
Ashtekar variables. Treating non-orthogonal boundaries was a matter of secondary concern. 
Thus, although certain elements of his discussion are similar to some of the developments 
of this paper, the focus is quite different. In particular he did not decompose the action J, 
and the decomposition of SI was with respect to variations of the Ashtekar variables rather 
than metric variables. He did not discuss background terms in detail and did not calculate 
any examples. 

In this paper we shall consider both a decomposition of Hayward's action and a decom- 
position of the variation of that action and show that they agree in their natural candidates 
for the quasilocal quantities. In doing this, we shall focus on the boundary lapse and shift 
functions rather than the full spacetime lapse and shift as was the case in JTTJ] . This will 
result in less complicated decompositions that also require less stringent restrictions on the 
comparison spacetime, in contrast with the approach in ref.Jll . We shall also argue that 
the boundary lapse and shift functions are the natural lapse and shift to consider. 

Before turning to those decompositions it is necessary to set out quite a few definitions. 
Those definitions will be the subject of section two. In section three we will perform the 
decompositions, examine the quasilocal quantities that naturally arise from those decom- 
positions, see how those quantities relate to conserved charges, and finally examine the 
background terms in some detail. Section four is made up of applications of the theory of 
section three. In that section we shall calculate the quasilocal energy (and other quanti- 
ties) seen by a spherically symmetric set of observers undergoing a variety of motions in 
Schwarzschild spacetime. 



2 Definitions 

Consider a region M. of an n-dimensional spacetime with metric tensor field g a p and on that 
region define a timelike vector field T a and a spacelike (n — l)-dimensional hypersurface 
So- This field and surface are sufficient to define a notion of time over A4. As a start, we 
let Eo be an "instant" in time. That is we choose to define all events happening on that 
surface as happening simultaneously. Next, consider a set of observers at locations G S 
(where the A index labels the individual observers). The past and future locations of these 
observers are uniquely determined if we specify that they must follow the paths x\{t) through 
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dx a 

spacetime where these paths satisfy the differential equation = T a subject to the initial 
condition x^(t ) = x^ . We then define "instants" of time to be t = constant surfaces. We 
label them E t and define the notion of past and future by saying that if t% < t 2 , then Si 
"happens" before E 2 . By construction the to hypersurface is E . Thus, from the vector field 
and original hypersurface we have imposed an observer dependent notion of time on our 
manifold according to the constructed time coordinate t. Next, we may break up T a into its 
components perpendicular and parallel to the E f by defining a lapse function iV and a shift 
vector field V a so that, 

T a = Nu a + V a , (1) 

where u a is defined so that at each point in Ai it is the future pointing unit normal vector 
to the appropriate hypersurface E t , and V a u a = 0. The lapse and shift then tell us how 
observers being swept along with the time flow T a move through space and time relative to 
the foliation. 

We specialize to the situation of interest for this paper. On the surface E define a surface 
Qo that is topologically an (n — 2) sphere. Qq locally bifurcates So - we will pick one of the 
regions as "inside" and the other as "outside" []. Now, propagate this surface through time 
according to the time flow T a . Then, by the local continuity of the time flow, fit will (at 
least locally) topologically remain an (n — 2) sphere in the hypersurface E t and still divide 
the "inside" from the "outside". Choosing time coordinates t\ and t 2 we define an (n — 1) 
dimensional timelike hypersurface B = {U t flt '■ ti < t < t 2 }. We then define M C M. as the 
region "inside" B bounded by the surfaces Eti and E i2 . E t and fl t represent foliations of M 
and B respectively. Figure [l] illustrates these concepts for a three dimensional Ai. 

We define unit normal vector fields for the various hypersurfaces. Already we have defined 
u a as the timelike unit normal vector field to the E t surfaces. Similarly, we may define u a as 
the forward pointing timelike unit normal vector field to the surfaces fl t i n the hypersurface 
B. The spacelike outward-pointing unit normal vector field to B is defined as n a . Then, 
by construction u a n a = and T a n a = 0. We further define n a as the vector field defined 
on B such that n a is the unit normal vector to fl t i n E 4 (fl t viewed as a surface in E t ). By 
construction u a h a = 0. 

We define the scalar field 77 = u a n a over B. If rj = everywhere, then the foliation 
surfaces are orthogonal to the boundary B (the case dealt with in refs. |3], 0), and the 
vector fields with the tildes are equal to their counterparts without tildes. We express u a 

1 Globally of course the "inside" and "outside" could be connected. Consider for example the case where 
So is a two torus, and f2o is a homotopically non-trivial circle in that surface. 

2 The definitions of u a and n a are consistent with ref. ^j, but interchanged with respect to those in ref. 

§■ 
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Figure 1: The Lorentzian region M, its assorted normal vector fields, and a typical element 
of the foliation. 



and n a in terms of u a and n a (or vice versa) as, 



n a = -n a - V u a and u a = \u a - rjn a , (2) 
A A 

or, 

fi a = \n a + r]u a and u a = \u a + i 1 n a , (3) 
A A 
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where A = t^s- 

Note too, that with T a n a = we may write, 

T a = Nu a + V a , (4) 

where we call N = \N the boundary lapse and V a = a^V 13 the boundary shift. 

Next consider the metrics induced on the hypersurfaces by the spacetime metric g a p. 
These may be written in terms of g a p and the normal vector fields. h a p = g a p + u a up 
is the metric induced on the S t surfaces, 7 aj g = g a p — n a np is the metric induced on B, 
and a a/ 3 = h a p — h a np = 7^ + u a xip is the metric induced on Q t . By raising one index of 
these metrics we obtain projection operators into the corresponding surfaces. These have the 
expected properties: h^u 13 = j^n 13 = a^n 13 = o a puP = 0, and h^h 13 ^ = h a y , 7°g7 /3 7 = 7 Q 7 , 
and a a p a^ = a a r 
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On choosing a coordinate system {x l ,x 2 ,x 3 } on the surface S we define h = det(h a p) 
(where in this case we take h a p as the coordinate representation of that metric tensor). 
We then map this coordinate system to each of the other E t surfaces using the time flow; 
combining this set of coordinates on each surface with the time coordinate x° = t we have a 
coordinate system over all of M. We define g = det(g a p). Similarly, choosing a coordinate 
system on Q t we define a = det(a" Qj g). Again, using the time flow to extend the coordinate 



system over all of B, we define 7 = det(j a p). It is then not hard to show |I]J that 



v /3^ = NVh and = Ny/a. (5) 

We also define the following extrinsic curvatures. Taking V a as the covariant derivative 
on M. compatible with g a p, the extrinsic curvature of E 4 in M. is K a p = — h J a h S gS/jUs = 
— \£ u hctp-, where £ u is the Lie derivative in the direction u a . The extrinsic curvature of B in 
M. is 9 Q /3 = — 7 7 Q 7 <S / 3V 7 n < 5, while the extrinsic curvature of fl t in is fc Q/ g = — cr 7 Q cr 5 /3 V 7 n5. 
Contracting each of these with the appropriate metric we define K = h af3 K a p, O = 7 a/3 Q/ 3, 
and k = a al3 k a /3. 

Finally, we define the following intrinsic quantities over M. and S t . On M, the Ricci 
tensor, Ricci scalar, and Einstein tensor are 7Z a/ 3, 7Z, and G a p respectively. On E t , D a is 
the covariant derivative compatible with h a p, while R a/ 3 and R are respectively the intrinsic 
Ricci tensor and scalar. 



3 Analysing the action 



For definiteness, we now take M. to be four dimensional. The generalization to other di- 
mensions is trivial. Given M C M. as described above and allowing for the inclusion of a 
cosmological constant, Hayward's action QTUJ is: 

d 4 x^/=q(n - 2A) + - / d 3 xVhK - - / d 3 xy/^Q 



1 

2k 



[ d A x^/^g{TZ - 2A) + - f d 3 xVhK-- [ 

J M K JT, k J E 

+ — [ d 2 x\fo sinh _1 (r]) — /, 
« 7n 



(6) 



where J s = J S2 _ 2l and J n = f n — J n , and if we choose a system of units where c—G—1, 
k = 8tt. I_ is a functional of the boundary metrics on dM. For simplicity, in the next two 



sections we will take 1 = 0, but in the section 3.4 we will allow it to be non-zero again 
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3.1 Variation of the action 



The variation of I with respect to the metric g a p is [ID 

51 = ±- [ d A xyf=7j(G a(3 + Ag aP )Sg a P (7) 

^ K J M 

d 3 x {P af3 5h aP ) + / d 3 x {ii a(3 Sj a(3 ) 
: J b 

+ / d 2 x f — sinh _1 (?7)5v / o : ) , 
Jn V K / 

where P aP = & (K afS - Kh a/3 ) and vr Q/3 = ^ (6 Q/3 - 9 7 a/3 ) . If we consider variations 
of the metric that leave the boundary metrics h a p and 7 Q /3 fixed, then all of the boundary 
terms are 0, and 81 = if and only if Einstein's equations hold over M. Thus / is the action 
that generates general relativity if we are considering variations of the metric over a bounded 
region of space such as M. 

Now, j a /3 is fully defined if we specify N, V a , and a a p. Thus variation of 7 a/ g is equivalent 
to a variation of these quantities, and we may rewrite the B term in the above with respect 
to such variations. During this calculation we repeatedly make use of the fact that 5w a ||ii a , 
5n a || n a . This is true because these one forms are defined by the requirement that u a v a = 
n a w a = for all vector fields v a G TH t , w a G TB. The metric does not figure in the 
definition of T£ t or TB nor is a metric required to calculate the action of a one form on 
a vector, so u a and n a are defined up to a normalizing factor independently of the metric. 
Thus, a^Sup = a aP 5u p = a a/3 5n p = a a0 8n p = 0. 

Expressing u a and h a in terms of u a and n a and writing 7 a/3 = <T a/3 — u a U/3, we have 

(e«/J _ 7 q / 3 )5 7q/3 = (6 Q/3 - Qo aP )5o a p + 2Q aP u a 5u? - 2QuJu a . (8) 

In the meantime, we may decompose Q a p into its parts that are perpendicular and parallel 
to the Q t to obtain 

©a/3 = k a/3 + 2u ia a f3 l ) u s V- y n s + UaU^rCa^), (9) 

where k a/3 = —a^aJV^ng (the extrinsic curvature of fl t in a hypersurface perpendicular 
to B) and a a = -u^V ' pu a is the acceleration of normal vector u a along its length. Also, on 
contracting Q a/3 with 7^, 

= k-n a a a , (10) 

where k = (T a ^k a p. Putting these three results together, a few lines of algebra 

produces 

(0«/3 _ &ff*)5 lafi = (k al3 - [k - n 7 5 7 ]a a/3 )<5a a/3 + 2a 7 n 5 W 1 u & 6u p - 2ku p 5u p . (11) 
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To complete this deconstruction, recall that the time-flow vector field is defined independent 
of the metric. Therefore ST a = which in turn implies that 6u a = — -~ (^u a 5N + 5V a ^j. 
Applying this to equation ([O]), substituting the result back into the variation of the action 
(0), and recalling fl5|) we obtain the final result 

SI = 0- [ d 4 x^(G a(3 + Ag af3 )6g af3 (12) 



-J d 3 x (P al3 5h a(3 )+ J d 2 x (^w&-\ri)8y/a\ 
— I dt I d 2 xy/a 

J JQt 



where e = -k, j a = — -a£n 5 V pus, and s a P = ^(k alS — [k 

From this result we can make a couple of useful observations. First, examining the 
initial and final hypersurfaces Si and E2 and their boundaries fl± and Q2 we see that P af3 is 
the Et hypersurface momentum conjugate to h a/3 while ^sinh -1 ^) is the f2 t hypersurface 
momentum conjugate to \fa. Secondly we see that —yfae is conjugate to the boundary 
lapse N, \fa] a is conjugate to the boundary shift V a , and \N^s af3 is conjugate to the 
boundary metric a a p. Following the Hamilton- Jacobi analysis of we identify e, j a , and 
s af3 as surface energy, momentum, and stress densities. If 77 = 0, these quantities coincide 
with those defined in || . Also note that each of these terms is explicitly independent of -q. 
They are defined with respect to the foliation of B only. 

3.2 Decomposing the action 

We now decompose / with respect to the foliation. To start, the 1Z — 2 A term of @ may be 
rewritten as 

/ d 4 xy/=g~(K - 2A) = [ d A x^g (R - 2A - K 2 + K af3 K a(S ) (13) 
Jm Jm 

-2 / d 3 x\fhK-2 / d 3 

X V~1 {Kr) + n a a a ) , 

Jt, Jb 

where a a = -u^V pu a is the acceleration of the foliation's unit normal vector field along its 
length. Next, if we substitute expressions for n a and u a from eq. @) into (|T0| ) then it is a 
simple matter to show that, 

Q + r]K + n a a a = \k + \u a V a r]. (14) 
A 
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Combining these two results the following expression for I is obtained. 

/ = _L f d A x^-g{R-2K~K 2 + K aP K^) (15) 

IK J M 

-- [ d 3 xy/^f ( ^ + \u a V a rj ] + - [ (fx^smh-^rj). 
kJb V a / K Jn 

Next, we apply the Einstein constraint equations. These are 



n = 



^G a pu a u? = -^5 (R - 2A + K 2 - K Q(3 K al3 ) = 0, and 



K 



H a = ^h£G^ = ^(D,Kf-D a K)=0. 

K K 



(16) 
(17) 



Combining these constraints with the Lie derivative definition of the extrinsic curvature 
K a/3 = —^£ u h a [3 = — (£rh a /3 — 2-D( a V^)) of E t in A4, we may rewrite the integrand of the 
remaining bulk term with respect to these constraints, a time derivative of the hypersurface 
metric, and a total divergence term: 



R - 2A - K 2 + K^JC* = ^=P aP £ T h aP ~H- -^=V°"H a -^D a 

Vh Vh \h9 N 



.Vh 



P af3 V B 



18) 



where P a/3 is the hypersurface momentum for S t that we discussed above. Then, using Stokes 
theorem on the hypersurfaces S t to move the total divergence term out to the boundaries 
fit and applying (||) we may write the action as: 



(19) 



/ = f d 4 x(P al3 £ T h al3 - NH-V a H a ) 

JM 

H — / d 2 X\fa sinh~ 1 (^). 
« J si 

Up to this point we have been working with the foliation of M and therefore with the lapse 
N, shift V a , and normal vectors u a and h a . On the term evaluated on B we now switch to 
work with the foliation of B and therefore the boundary lapse N, the boundary shift V a , 
and normal vectors u a and n a . Then, 



ATA; = ^-Nk-7]Na aP V a ui3 
A 



(20) 



-h a (K af3 - Kh aP ) V p = Nrja a ^V a U(3- Nrfl + n a V^W p u a + AV^V^, (21) 



9 



where k = a al3 k a p. Combining these two results with (ffl), which can be used to show that 

/ (fx^asmhT 1 ^) - dt d 2 xV^XNu a V a r] (22) 
Jn J Jn t 



= J dt J d 2 x ((£ T y/a) sinh- 1 (ri) + ^/aXV a V a rj^ 

we obtain the following decomposition of the action, 

/ = f d 4 x (P al3 £ T h a ^ -NH- V a H a ) + - / dt [ d 2 x{£ T ^a) sinh -1 ^) (23) 
Jm K J Jn t 

- J dt J d 2 x^ (fie - V a jS 



where e and j a are the energy and surface momentum densities that we obtained from the 
variational calculation. 

The terms of this expression will be familiar to anyone who is familiar with refs. || ^ [TTJ. 
Specifically, e and j a are exactly the energy surface density and momentum surface density 
that the observers on the boundary would measure if the foliation of M were perpendicular 
to B. A little thought shows that these quantities are the ones that would be reasonable for 
observers restricted to surface B to measure. Such observers are cognizant of the foliation 
of the boundary (for the foliation has been defined to correspond to their notion of simul- 
taneity), but being restricted to the surface they have no way of associating that foliation 
with a foliation of M as a whole. Viewed another way, there are no observers in the interior 
of M and therefore no unique way to extend the "instants" of time into that region. As 
such, it does not seem to physically make sense for the observers to measure the energy and 
momentum surface densities with respect to the foliation S 4 , the lapse N, and the shift V a 
that we have defined but they cannot observe. Rather, as observers travelling along B they 
would naturally (locally) extend the foliation of B into a foliation that is perpendicular to 
T a - and therefore in effect be considering a foliation S t (locally) defined around B with 
lapse iV and shift V a . Then they would measure the quantities that we have found naturally 
arise from the action. 

We may also define a Hamiltonian. In elementary classical mechanics with one degree of 
freedom, the action I and Hamiltonian H are related by the equation I = pq — H, where q is 
the variable giving the configuration of the system and p = |^ is the conjugate momentum. 
Extending this to the system under consideration |J, h a p and yfa are configuration variables 
while P a/3 and ^sinh _1 (?7) are their conjugate momenta and so the Hamiltonian is: 

H= [ d 3 x{NH + V a H a ]+ [ d 2 x^{Ne~ V a j a ). (24) 
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Again this quantity is indifferent to the intersection angle between the foliation of M and 
the boundary. For solutions to the Einstein equations, it is defined entirely with respect to 
the foliation of B. Note that this Hamiltonian does not agree with that proposed in |TT] 



where the problem was approached from the point of view of the foliation of M rather than 
that of B. 



3.3 Conserved Charges 

The discussion of conserved charges presented in H carries over exactly into this work. Thus 
if £ a e TB is a vector field in the boundary B and £^ja/3 = (ie. it is a Killing vector field), 
then we may define an associated conserved charge 

Q € = ( d 2 x^i a {eu a +j a ). (25) 
Jn 

If T a is a Killing vector field then the Hamiltonian H as defined above is a conserved charge. 
If there is an angular Killing vector field <p a G Tfl then the angular momentum 

J+= [ d 2 x^ a J a , (26) 
Jn 

is also a conserved charge. 

3.4 Background Terms 

We now return to the reference term J. Defined as it is as a functional of the boundary 
metrics, it is clear that for a metric variation that leaves the boundary metrics unchanged, 
5I_ = - therefore its exact form does not affect the equations of motion. This degree of 
freedom in the definition of / may equivalently be viewed as the freedom to define zero points 
of the energy, momentum, and Hamiltonian. Specifically, it allows us to choose a reference 
spacetime for which we wish these quantities to be zero. For asymptotically flat spacetimes 
we would normally choose Minkowski space as the reference spacetime, but other choices 
may be made if we are studying spacetimes with other asymptotic behaviours - for example 
asymptotically anti-deSitter space ||. 

Given a reference spacetime (M,q g ), we embed (ft, cr a p) in that spacetime and define a 
vector field T a over the embedded (T2, a^p) such that T a T a = T a T a and the components of 
£ro_ al3 = £T&af3 0- These conditions ensure that the boundary lapse and the components of 

3 We leave aside the issue as to whether this is possible in all cases. We will consider several examples 
where it is but in general an arbitrary surface cannot be embedded in an arbitrary higher dimensional space. 
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the boundary shift vector as calculated from T_ a are equal to those calculated for T a in the 
original spacetime. We then define 

dt [ d 2 x^fa[Ns-V a j ], (27) 

where e and J q are defined in the same way as before except that this time they are evaluated 
for the surface Q embedded in the reference spacetime. Thus, the net effect of including / is 
to change e — > e — e and j a — > j a — j . 

Physically these conditions correspond to demanding that an observer living in the surface 
Q and observing only quantities intrinsic to that surface (as it evolves through time) cannot 
tell whether she is living in the original spacetime or in the reference spacetime. From 
another point of view the observers have calibrated their instruments so that they will 
always measure the quasilocal quantities to be zero in the reference spacetime - no matter 
what kind of motion they undergo. 

This definition of / differs slightly from both the one used in || and the one used in 
TTJ. In the former case Q was embedded in a reference three dimensional space and no 



demand was made of £xO'ai3', however in all examples considered in ref. || (and indeed 
in the subsequent work) £T&ap = 0, so insofar as that formalism has been pursued within 
the literature it agrees with the formalism considered here. Note that if we do not include 
the conditions on £TO~ a p then boosted observers in Minkowski space will observer non-zero 
quasilocal energies which is clearly an undesirable situation. 

In ref. [[TI| (S,7 a/ g) as a whole is embedded in the reference four space (M,g>). That 
requirement is essentially the global version of our definition of / and as such will locally yield 
the same results as our definition though it is somewhat harder to apply computationally. 
Beyond that condition they further require that the reference spacetime be foliated in such 
a way that rj in the reference spacetime is the same as in the original spacetime. Such a 
condition is neither necessary nor desirable in our approach which doesn't concern itself with 
the foliation of the spacetime as a whole. Finally we note that in the approach used in ref. 
[ ]"T| the inclusion of this background term is necessary to remove an 7] dependence in the 



Hamiltonian - this dependence does not occur in our approach. 



4 Examples 

We now consider some sample calculations. For simplicity we will work with static spherically 
symmetric spacetimes parameterized with the natural spherical coordinates {t, r, 8, </>} and 
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therefore with metric 



ds 2 = -F(r)dt 2 + -^j-y + r 2 (d9 2 + sin 2 9d(j ) 2 ). 



(28) 



In each case we will consider a surface of observers fl denned by r = r and t = t . Then 
geometrically fl it is a two sphere with metric 

ds 2 = r 2 (d9 2 + sin 2 9d(f) 2 ), (29) 

(where we have parameterized fl with the same 9 and cf> as the full space). If we then consider 
the timelike unit vector field, 



N 



R 2 

i + ^ — ,i?,e,$ 

F N 2 F 



(30) 



where R = i?(r, 9, </>), © = 6(r, 0), and $ = $(r, 0) are general functions of r, 6>, and </>, 
the boundary lapse and boundary shift functions are easily found to be 

N 2 = 1 + r 2 (9 2 + sin 2 #$ 2 ), 

and 

V a = [0,0,6,$], 

while 



u a = 



1 R 2 R 

ll + -= — ,^,0,0 
FV N 2 F N 



and 



Then, a straightforward calculation yields 



2 L R 2 

e = \ F + , and 

«T V N 2 



Ja 



0,0, 



9 {R\ d (R 



89 \ N J dtp V N 



Finally, we calculate the Q components of £r a 'ai3 as 

(£t^^)^ = 2rsin# (i?sin# + OrcosO + rsin 



2r { R + r m)' 



<9$ 



(31) 
(32) 

(33) 



(34) 

(35) 
(36) 



(37) 
(38) 
(39) 
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Note that there is no F dependence in any of these components. Thus if we wish to calculate 
quasilocal quantities for observers moving through Schwarzschild space using Minkowksi 
space as a reference spacetime (as in the following examples), on embedding Q (which is 
trivial for a sphere) and setting 



(40) 

the metric a a p and its derivative £T®otp will be the same for both Schwarzschild and 
Minkowski space. 

We now specialize to specific examples using the system of units where k = 8n. 



R 2 

T a = N\ll + ^,R,e,® 
N 2 



4.1 Static Observers 

For our first example, we'll consider a spherical set observers holding themselves static with 
respect to a Schwarzschild black hole (F — 1 — — ) and take flat Minkowski space (F = 1) 
as our reference spacetime. Then for both spacetimes i? = O = $ = 0, N = 1, and V a = 0. 
Substituting these data into the above expressions we obtain 




Then, the total measured energy is 

E = j d 2 x^{e - 1) = r (l - yjl-^j . (42) 

This is the standard result as obtained in || . Taking the limit as r — > oo we obtain E — > m 
as would be expected, while asr-» 2m (the Schwarzschild horizon), E — > 2m. 

With N = 1 and the shift vector 0, the Hamiltonian H = E. Since T a is a Killing 
vector in this case, H is a conserved charge. There is also a conserved angular momentum 
associated with each of the regular three spherical Killing vectors. j a is zero however, so 
each of these charges vanishes. 

4.2 Radially Infalling Observers 

A more interesting example is the case of observers taking their measurements as they fall 
radially along geodesies towards a Schwarzschild hole. Such motion is described by solutions 
to the geodesic equation. For observers who were stationary as they started falling in from 
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infinity, the geodesic equation reduces to ^ = — \/^r, where r is the proper time coordinate. 



Then, R = -y^, G = $ = 0, N = 1, again \/ a = 0, and 
The total measured energy is, 

/ / 9m \ 

(44) 



E = j d 2 x^a(e - e) = r {^j 




As for static observers as r — > oo, £" — > m. Of course, this is not really surprising since the 
radially infalling observers at infinity actually are static! Over the rest of the range the two 
energy measures are not the same. In particular as r — > 2m, E -> 2m( v / 2 - 1). 

As for the first example, the momentum terms are zero and the N — 1, and so H = E. 
T a is no longer a Killing vector however, so this is no longer a conserved charge. Physically 
of course this is to be expected since the observers are moving radially inwards and therefore 
through the gravitational field. As time passes therefore the amount of gravitational field 
energy contained within Q t changes. Again the three angular momenta are conserved but 
each has the uninteresting value of zero. 

4.3 Radially Boosted Observers 

We next consider a set of observers who are boosted to travel radially with "constant" 
velocity v. By constant velocity here we mean that a second set of observers dwelling on a 
t = constant surface and being evolved by the timelike vector field [1,0,0,0] will measure 
the first set as hav ing velo city v and acceleration 0. 

Then R = ^Jl — ^ (7 = - the standard Lorentz factor from special relativity), 

O = $ = 0, N — 1, and once more V a = 0. A simple calculation then obtains, 




^-i=4Z7"! \ ' 1 \ ' —I (15) 

and the total measured energy is, 



E= f d 2 xy/a(e - e) = 7 r ( Jl - — - Jl - — | . 1 16) 

Jn 



r v r 



In this case as r — > 00, E — > y \ — v 2 m, while as r — > 2m, £" — > 2m. 
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Once more with N = 1 and V a = 0, H = E. Again T a is not a Killing vector, so this is 
not a conserved charge. The angular momenta are conserved charges though again each is 
zero. 

At first glance it may seem unusual that in the r — > oo limit E oc — . Extrapolating 
from special relativity we would perhaps expect E oc 7m. Physically however, it is clear 
that there is a flow of gravitational field energy through the surface Q. That is, there is a j\- 
component of the momentum. This momentum may be seen as "drawing off" some of the 
energy. We will not investigate the issue further in this paper, though it is addressed in the 



last example of [^T| to a certain extent by the invariant quantities defined in |L2 



4.4 Z-Boosted Observers 

Finally we consider a set of observers who are boosted to travel "in the ^-direction" with 
"constant" velocity v. By constant we again mean with respect to other observers who are 
dwelling on t = constant surfaces and being evolved by the timelike vector field [1,0,0,0]. 



Then, R = jvcos9Jl-^, = ^f^, $ = 0, N = y/l + 7 V sin 2 9, and V c 
[0,0,2^,0]. We now have 



1 1 L 2mv 2 cos 2 6 . , 

^=7-^^^ \/i : \/!-— . (47) 



47rr Vl - v 2 cos 2 9 

and while = 3 = 0, 




v sin 9 I 1 — 



~ 3e & 8tt(1 - v 2 cos 2 9) 1 y 1 - w } os2g 1 I ' 

In this case, the e — e doesn't integrate into a nice tidy form as it did in previous examples. 
Instead we will consider the two usual limiting cases. For r — > 00, 

(e - e) r _ = 008^(9 and ( Jq - j Q j_ = . (49) 

Then, integrating over the {t = 0, r — >• 00} two surface we obtain, 

Jn 



e - e 



^(7T^+^),aod (50) 



i/^ = J d 2 x^(N(l-l)-V a (3 a -3j^ 



vl — f 2 m. (51) 
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Meanwhile for r — > 2m, 

1 /_ ~ \ vsin6 

(£ ~ " = AM2m) ^ (*■ - ^L,™ = 4.(1 -ttWfl) ' (52) 

Integrating over the {t = 0,r = 2m} two surface we obtain, 

E 2m = / d 2 xy/a(i - |) 
in 

= 2m, and (53) 
H 2m = J d 2 x^(N(e-e)-V a (j a -jj' ) ) 

( 7 2 \ 

= m 1 — 27 H — arcsinf H arctanh v . (54) 

\ v l v ) 

As in the previous cases T a is not a Killing vector and so H is not a conserved charge, 
as we would physically expect. Note that in this situation we have a non-zero component 
of j a . Despite this we still do not have any non-zero conserved angular momenta. To see 
this recall that the three linearly independent spherical Killing vectors are 0° = [0,0,0, 1], 
02 — [0, 0, sin 0, cos cot 9], and 03 = [0, 0, cos0, — sin 0cot 9}. Then 4>\3a = 0, (fi^ja is 
proportional to sin 0, and 0f is proportional to cos 0. Therefore, as would be expected 
physically and by symmetry, on integration = J^ 2 = J^ 3 = 0. 

The comments made in the previous example regarding the fact that = ^ apply here 
as well. 

5 Discussion 

In this paper we have seen that the orthogonality restriction of || may be lifted without too 
much difficulty Further, by concentrating on the foliation of the boundary B rather than 
the spacetime region M we avoid many of the technical complications of the non-orthogonal 



treatment of [|11[] , and obtain definitions of quasilocal quantities that are manifestly indepen- 
dent of the intersection angle between the foliation of M and the boundary B independent 
of our choice of the background spacetime. 

In our choice of how to calculate the reference term / on the background spacetime we 
have given local conditions that modify those of || in a way that is more appropriate if we 
are considering moving observers. These conditions at the same time remain simpler and 
easier to implement than those required in ||11|| . 

In general the Hamiltonian and quasilocal quantities such as angular momenta are de- 
pendent on the motion of the observers as we have seen in several examples. One somewhat 
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counterintuitive observation is that the observed mass of a source decreases rather than in- 
creases with the motion of the observers who are measuring that mass. This is a consequence 
of choosing our observers in such a way that there is a net flow of gravitational field energy 
through the surface Q. 
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